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IS  ABSTRACT 


■Experiments  were  conducted  in  the  northern  Gulf  of  Mexiro  (30°00'06"N,  85°54'03"W) during 
July  1969  and  March  1970  to  study  propagation  characteristics  of  acoustic  normal  modes.  This  shallow- 
water  area  had  a  constant  depth  (31  i  2.5  m)  over  a  propagation  path  of  at  least  20  km  along  a  bearing 
line  of  285°T.  Short-pulse,  3  cycle  signals  were  emitted  from  400-  or  750-Hz  center  frequency  sources. 
The  vertical  distributions  of  pressure  amplitude  of  each  of  the  two  lowest  order  normal  modes  of  these 
signals,  measured  by  a  vertical  line  of  hydrophones,  agreed  well  with  distributions  calculated  by  computer 
using  a  two-layer  fluid  model.  Temporal  variations  of  modal  intensity  averaged  over  depth  for  both  the 
first  and  the  second  modes  were  observed  for  7-min  runs. ,  Although  the  specific  causes  of  these  variations 
were  not  determined,  there  were  indications  that  the  vertical  distributions  of  pressure  amplitude  remained 
comparatively  stable  and  signal  fluctuations  were  occurring  uniformly  with  depth.  Time-averaged  attenua¬ 
tion  coefficients  were  also  obtained  for  the  first  and  second  modes  for  7-mln  runs.  When  a  mode  was 
predicted  to  interact  strongly  with  the  bottom,  its  measured  attenuation  was  high  compared  to  cases  where 
that  mode  was  predicted  to  interact  weakly  with  the  bottom. 
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ABSTRACT 

Experiments  were  conducted  in  the  northern  Gulf  of  Mexico  (30°  00'  06" 
N,  85°54'03"W)  during  July  1969  and  March  1970  to  study  propagation 
characteristics  of  acoustic  normal  modes.  This  shallow-water  area  had  a 
constant  depth  (31  ±  2.5  m)  over  a  propagation  path  of  at  least  20  km  along 
a  bearing  line  of  285°T.  Short-pulse,  3  cycle  signals  were  emitted  from  400- 
or  750-Hz  center  frequency  sources.  The  vertical  distributions  of  pressure 
amplitude  of  each  of  the  two  lowest  order  normal  modes  of  these  signals, 
measured  by  a  vertical  line  of  hydrophones,  agreed  well  with  distributions 
calculated  by  computer  using  a  two-layer  fluid  model.  Temporal  variations 
of  modal  intensity  averaged  over  depth  for  both  the  first  and  the  second 
modes  were  observed  for  7-min  runs.  Although  the  specific  causes  of  these 
variations  were  not  determined,  there  were  indications  that  the  vertical  dis¬ 
tributions  of  pressure  amplitude  remained  comparatively  stable  and  signal 
fluctuations  were  occurring  uniformly  with  depth.  Time-averaged  attenuation 
coefficients  were  also  obtained  for  the  first  and  second  modes  for  7-min 
runs.  When  a  mode  was  predicted  to  interact  strongly  with  the  bottom,  its 
measured  attenuation  was  high  compared  to  cases  where  that  mode  was  pre¬ 
dicted  to  interact  weakly  with  the  bottom. 
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SOME  MEASUREMENTS  OF  ACOUSTIC  NORMAL  MOPS 
PROPAGATION  IN  SHALLOW  WATER 


INTRODUCTION 

Validation  tests  of  normal  mode  acoustic  propagation  theory  have  been  conducted  in 
a  shallow-water  area  of  the  Gulf  of  Mexico.  Acoustic  sources  were  mounted  on  a  fixed 
platform,  a  facility  of  the  Naval  Ship  Research  and  Development  Center,  located  approx¬ 
imately  20  km  offshore  from  Panama  City,  Florida,  in  a  water  depth  of  31  m.  The 
acoustic  field  was  sampled  by  a  vertical  line  of  hydrophones  suspended  from  the  USNS 
GIBBS  (T-AGOR-1). 

The  objective  was  to  compare  the  vertical  distribution  of  pressure  amplitude  in  indi¬ 
vidual  normal  modes  of  propagation  with  those  predicted  by  an  idealized  wave  theory 
model  (1)  of  the  area.  The  field  of  short,  gated  sine  wave  pulses  was  sampled  from  near 
the  surface  to  near  the  bottom  at  several  ranges  along  a  fixed  bearing  line.  Since  the 
group  velocities  of  signals  propagating  in  discrete  normal  modes  are  in  general  different 
from  mode  to  mode,  the  modal  arrivals  in  some  cases  could  be  separated  by  time  of  ar¬ 
rival.  In  those  cases  where  individual  modes  were  resolvable,  the  vertical  distributions  of 
pressure  amplitude  were  measured  and  compared  with  the  theory,  in  addition,  measure¬ 
ment  of  the  modal  field  intensities  permitted  the  calculation  of  mode  attentuation 
coefficients. 

Two  experiments  were  conducted  over  the  same  propagation  path;  one  was  in  July 
1969  and  the  other  in  March  1970.  In  the  first  experiment  a  frequency  of  400  Hz  was 
used  and  the  vertical  line  of  hydrophones  contained  eight  uniformly  spaced  elements. 

The  prevailing  sound  speed  profile  bad  a  negative  gradient.  During  the  second  experiment 
a  frequency  of  750  Hz  was  used  and  the  number  of  elements  in  the  line  of  hydrophones 
was  increased  to  twelve.  The  variety  of  observed  sound  speed  profiles  included  positive, 
zero,  and  negative  gradients. 

A  previous  report  (2)  described  the  analytical  model  and  the  experimental  technique 
used  in  these  experiments  and  presented  the  mode  attenuation  data  obtained  during  the 
first  experimental  period.  The  present  report  documents  the  results  obtained  from  the 
second  set  of  experiments  together  with  some  additional  analysis  of  the  July  1969  data. 
Included  are  calculations  of  attenuatior  coefficients  for  the  first  and  second  modes  at 
750  Hz  and  for  the  first  mode  at  400  Hz,  and  a  comparison  of  measured  and  calculated 
pressure  amplitude  distributions  for  these  two  modes  at  both  400  and  750  Hz.  A  spatial 
filtering  technique,  based  on  the  orthogonal  property  of  normal  modes,  is  described.  The 
output  of  such  a  filter  is  used  to  calculate  the  modal  attenuation  coefficients  and  to  ex¬ 
amine  the  temporal  stability  of  the  modal  intensities  averaged  over  depth. 
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ANALYTICAL  MODEL 
Theoretical  Formulation 

The  model  used  to  predict  the  modal  vertical  pressure  amplitude  distribution  and 
group  velocities  assumes  the  following  environmental  constraints:  (a)  the  water  layer 
depth  H  is  constant  along  the  propagation  path;  (b)  the  water  layer  has  a  constant  density 
Pl  and  a  sound  speed  c1  ( z )  which  is  only  a  function  of  the  depth  z;  and  (c)  the  bottom 
is  a  semi-infinite  fluid  of  constant  density  p2  constant  sound  speed  c2.  The  geometry 
of  this  model  is  shown  in  Fig.  1. 
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Fig.  1  —  Geometry  of  two-layer  analytical  model 


Tire  pressure  amplitude  Pn  of  the  nth  mode  established  by  a  pulsed  signal  emitted 
from  an  omnidirectional  point  source  of  center  angular  frequency  to  and  unit  source 
strength  positioned  at  a  depth  z0  can  be  related  to  the  free-field  source  level  S,  expressed 
in  decibels  relative  to  1  dyne/cm2  at  1  m  by  (2), 


20  log  Pn  (r,  £)  =  S+20  log 


+  10  log 


Here 

n  is  the  mode  number 
r  is  the  range  from  the  source 
f  =  z/H  is  the  normalized  depth 
kn  is  the  horizontal  component  of  the  wave  number  for  the  nth  mode 
T0  is  the  pulse  length  at  the  source 
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ATn  is  the  change  in  pulse  length  for  the  nth  mode  at  range  r  (due  to  dispersion) 
an  is  the  attenuation  coefficient  of  the  nth  mode. 

un  (?o)  anc*  un  (f)  are  the  eigenfunctions  of  the  nth  mode  evaluated  at  the  depth  of 
the  source  and  the  depth  of  the  receiver,  respectively. 

The  first  two  terms  on  the  right  side  of  Eq.  (1)  represent  the  equivalent  source  level  for 
the  nth  mode.  The  second  term,  the  mode  excitation  function,  expresses  the  decrease  in 
the  nth  mode  signal  level  from  the  free-field  level  S.  A  typical  value  of  this  decrease  for 
these  experiments  is  —20  dB.  The  third  term,  10  log  r,  represents  cylindrical  spreading. 
The  fourth  term,  10  log  T0/(T0  +  ATn),  is  an  approximation  of  the  decrease  in  peak  pres- 
sure  amplitude  in  a  short  pulse  caused  by  the  dispersive  lengthening  of  the  pulse.  The 
last  term,  anr,  represents  the  modal  attenuation  at  range  r.  Attenuation  results  from  a 
number  of  mechanisms  such  as  absorption  in  the  water,  scattering  by  the  boundaries  and 
by  inhomogeneities  in  the  water,  generation  of  shear  waves  in  the  bottom,  and  absorption 
in  the  material  of  the  bottom.  In  the  present  stage  of  development  of  the  analytic  model, 
only  the  summed  effect  of  all  loss  mechanisms  associated  with  interaction  with  the  bottom 
is  treated. 

Of  the  factors  occurring  in  Eq.  (1),  kn  and  un  (J)  were  calculated  by  a  normal  mode 
computer  program  using  an  iterative  finite  difference  technique  (3).  The  input  parameters 
of  this  program  were  the  water  and  bottom  densities,  and  p2,  respectively;  the  water 
sound  speed  profile  Cj  (z);  the  bottom  sound  speed  c2;  the  water  depth  H;  and  the  source 
center  angular  frequency  w.  This  program  also  calculated  the  mode  group  velocities  Un, 
given  by 


U„ 


dcj 


Influence  of  Sound  Speed  Profile  on  Modal  Propagation  Characteristics 


(2) 


A  variety  of  vertical  sound  speed  profiles  was  observed  during  the  March  1970  ex¬ 
periment.  Profiles  a,  c,  and  j  (see  Tables  1  and  2)  are  measured  examples  of  negative, 
zero,  and  positive  gradient  profiles,  respectively.  Their  influence  on  modal  propagation 
characteristics  maybe  examined  for  a  fixed  source  center  frequency  (in  this  case,  750  Hz). 


Vertical  pressure  amplitude  distributions  calculated  for  the  first  four  modes  using 
the  zero  gradient  profile  (Fig.  2b)  suggest  that  the  modal  pressures  are  approximately 
symmetrically  (odd-numbered  modes)  or  antisymmetrically  (even-numbered  modes)  dis¬ 
tributed  about  the  center  of  the  water  column.  The  distributions  associated  with  the 
negative  and  positive  gradients  (Figs.  2a  and  2c,  respectively)  illustrate  that  the  first  mode 
(and  only  the  first  mode)  has  a  majority  of  its  energy  confined  to  propagate  in  the  slower 
sound  speed  water  layer.  A  mode  is  considered  “confined”  to  a  slower  sound  speed 
water  layer  when  its  calculated  pressure  amplitude  function  has  an  inflection  point  at  a 
depth  where  that  function  has  a  nonzero  value.  Such  a  mode’s  wave  number  kn  will  be 
greater  than  or  equal  to  the  minimum  value  of  the  total  wave  number  K  possible  in  the 
water  column  (i.e.,  K  -  w/maximum  value  of  cj  ({’)).  In  general,  constriction  of  the 
modal  pressure  amplitude  distribution  requires  sharper  gradients  to  confine  more  and 
higher  order  modes. 
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Table  1 


Measured  Sound  Speeds  (in  m/sec)  for  March  1970  Experiment 


Depth 

ra 

1 - - - -  - - 

Velocimeter  Cast 

a 

b 

c 

d 

e 

£ 

! 

g 

mm 

i 

j 

0 

1513.9 

1512.9 

1508.7 

1508.5 

1509.4 

1507.0 

1508.3 

1507.0 

1508.2 

1505.8 

2 

1513.9 

1512.9 

1508.7 

1508.5 

1507.0 

1508.3 

1507.0 

1505.8 

4 

1513.9 

1512.9 

1508.8 

1508.5 

1509.4 

1507.0 

1508.4 

1507.1 

1506.4 

1505.8 

6 

1513.9 

1512.9 

1508.8 

1508.6 

1509.5 

1507.3 

1508.5 

1507.2 

1506.4 

1505.9 

8 

1514.0 

1513.0 

1508.8 

1508.6 

1509.6 

1507.3 

1508.4 

1507.2 

1506.7 

1505.9 

10 

1514.1 

1513.1 

1508.8 

1508.7 

1510.1 

1508.1 

1508.9 

1507.5 

1506.9 

1506.1 

12 

1514.0 

1513.1 

1508.9 

1508.8 

1510.4 

1508.6 

1510.1 

1507.7 

1508.8 

1511.6 

14 

1513.8 

1513.3 

1508.9 

1508.8 

1511.1 

1509.0 

1511.5 

1508.7 

1511.9 

1512.2 

16 

1513.9 

1513.6 

1508.9 

1508.7 

1513.3 

1510.2 

1512.9 

1511.0 

1512.3 

1512.2 

18 

1512.0 

1513.4 

1509.0 

1508.8 

1513.5 

1510.2 

1513.1 

1512.8 

1512.4 

1512.3 

20 

1511.3 

1511.1 

1509.0 

1508.9 

1513.6 

1510.6 

1513.1 

1512.9 

1512.4 

1512.3 

22 

1511.0 

1510.1 

1509.3 

1509.0 

1513.6 

1511.2 

1513.1 

1512.9 

1512.5 

1512.4 

24 

1510.9 

1509.9 

1509.3 

1509.6 

1513.6 

1511.2 

1513.0 

1512.9 

1512.5 

1512.4 

26 

1510.9 

1509.9 

1509.5 

1509.5 

1513.6 

1511.3 

1513.0 

1513.0 

1512.5 

1512.4 

28 

1510.9 

1509.9 

1509.8 

1509.7 

1513.7 

1511.3 

1513.2 

1513.1 

1512.6 

1512.5 

30 

1510.9 

1510.0 

1510.0 

1510.4 

1513.7 

1511.4 

1513.2 

1513.1 

1512.6 

1512.5 

bottom* 

1510.9 

1510.0 

1510.0 

1510.4 

1513.7 

1511.4 

1513.2 

1513.1 

1512.6 

1512.5 

*The  sound  speed  at  the  bottom  of  the  water  column  was  estimated  by  extrapolation  of  measured  values. 
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Table  2 

Environmental  Data  Associated  with 
Sound  Speed  Measurements  in  Table  1 


Veloci- 

meter 

Cast 

March  1970 
Date-Time 
Group 

Range* 

(km) 

Wind 

Speed 

(knots) 

Wind 

Direction 

(T) 

Wave 

Height 

(m) 

a 

130000Z 

11.9 

23 

270 

KOI 

b 

121438Z 

7.9 

t 

t 

mm 

c 

161659Z 

9.9 

10 

075 

1.0 

d 

161427Z 

9.9 

20 

075 

1.0 

e 

131730Z 

7.9 

t 

M 

f 

170003Z 

19.8 

12 

160 

0.5 

g 

131952Z 

9.8 

26 

295 

mrm 

h 

132354Z 

13.7 

■s 

320 

i 

160354Z 

9.9 

Si 

350 

n 

_ L _ 

160033Z 

13.8 

18 

280 

1.0 

♦Range  along  a  285°T  bearing  line  from  Stage  I  (30°00,06,'N, 

85°54'03"W). 

fData  not  taken. 
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MOOE  I  MODE  2  MODE  3  MODE  4 

Fig.  2  —  Vertical  distribution*  of  pressure  amplitude  of 
the  first  four  mode*  at  750  Hz  were  calculated  using 
three  different  sound  speed  profiles  from  Table  1 : 
(a)  negative  gradient  profil  (b)  isovelocity  profile  c. 
and  (c)  positive  gradient  pro  j. 
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Group  velocity  calculations  (Fig.  3b)  using  the  zero  gradient  profile  indicate  that,  for 
the  frequencies  shown,  mode  group  velocities  decrease  with  increasing  mode  number.  Suf* 
ficiently  sharp  gradients  may  change  this  order  (see  Figs.  3a  and  3c)  by  confining  a  majority 
of  the  energy  of  otherwise  faster  low-order  modes  to  a  slower  sound  speed  layer. 


MODE 

NUMBER 


650  750  850 

FREQUENCY  (Hz) 

Fig.  3  —  Group  velocities  of  the  first  four  modes  were  calculated 
using  three  different  sound  speed  profiles  from  Table  1:  (a)  neg¬ 
ative  gradient  profile  a,  (b)  isovelocity  profile  c,  and  (c)  positive 
gradient  profile  j. 
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Since,  for  a  given  bottom  and  fixed  mode  number,  the  mode’s  pressure  amplitude  at 
the  bottom  increases  as  the  percentage  of  the  mode’s  total  energy  that  propagates  in  the 
bottom  increases,  it  seems  reasonable  to  assume  that  the  bottom  attenuation  of  a  mode 
(that  portion  of  its  total  attenuation  due  to  interaction  with  the  bottom)  would  increase 
as  the  mode’s  pressure  amplitude  at  the  bottom  increases.  For  a  zero*  or  positive-gradient 
profile  the  calculated  modal  bottom  pressure  amplitude  increases  with  increasing  mode 
number  (Figs.  2b  and  2c).  However,  for  negative  gradient  profiles  (Fig.  2a)  sufficiently 
sharp  to  confine  some,  but  not  all,  of  the  modes  near  the  bottom,  some  higher  order  un¬ 
confined  modes  may  liave  lower  bottom  pressure  amplitudes  than  lower  order  confined 
modes.  The  bottom  pressure  amplitude  of  the  first  mode  is  highest  for  a  confining  nega¬ 
tive-gradient  profile  (Fig.  2a),  intermediate  for  an  unconfining  zero-gradient  profile  (Fig. 
2b),  and  lowest  for  a  confining  positive-gradient  profile  (Fig.  2c).  Such  trends  have  been 
observed  for  higher  modes  when  sufficiently  strong  profiles  to  confine  them  were  used. 
Qualitative  evidence  that  the  measured  total  mode  attenuation  follows  the  expected  atten¬ 
uation  trends  for  the  first  and  second  modes  associated  with  profiles  a,  c,  and  j  will  be 
presented  in  the  Mode  Attenuation  Coefficients  section. 


EXPERIMENTAL  PROCEDURES 

This  section  describes  details  of  the  March  1970  experiment.  Similar  information 
pertaining  to  the  July  1969  experiment  is  found  in  Ref.  2.  The  bottom  profile  along  a 
285°  true  bearing  line  from  Stage  I  (30° 00' 06"  N,  85°54'03"W)  was  measured  with  a 
shallow-water  fathometer  over  ranges  from  0.9  to  17.8  km  from  the  Stage.  (Figure  4 
depicts  the  geography.)  Readings  indicated  an  average  depth  of  31.3  m  with  maximum 
bottom  contour  variations  of  ±  2.5  m.  Tidal  depth  changes  measured  at  Stage  I  were 
approximately  0.7  m.  The  source  was  rigidly  mounted  on  a  carriage  that  could  be  lowered 
and  raised  on  rails  attached  to  Stage  I.  A  vertical  line  of  12  hydrophones  spaced  at  2-m 
increments  was  lowered  over  the  side  of  a  ship  to  depths  from  either  5  to  27  or  6  to  28 
m.  Each  hydrophone  was  equipped  with  a  preamplifier  having  16.3-dB  gain.  The  mechan¬ 
ical  Q  of  the  750-Hz  omnidirectional  transducer  was  approximately  4,  so  acoustic  signals 
did  not  build  up  to  a  maximum  possible  output  level  with  the  three-cycle  electrical  input 
pulse.  The  maximum  buildup  of  the  source  output  corresponded  to  a  source  level  of 
103.2  dB  relative  to  1  dyne/cm2  at  1  m. 

Individual  runs  lasted  about  7  min  and  consisted  of  approximately  60  three-cycle 
pulses.  Data  were  recorded  and  stored  by  an  FM  analog  tape  recorder  and  subsequently 
rerecorded  by  an  oscillographic  recorder  for  processing.  A  block  diagram  of  the  recording 
instrumentation  is  shown  in  Fig.  5. 

Sound  speed  profiles  were  recorded  two  or  three  times  daily  using  two  adjacently 
mounted  velocimeters.  They  were  lowered  to  the  bottom  to  dislodge  any  trapped  air 
bubbles  and  then  raised  in  2-m  increments  to  the  surface.  Table  1  contains  all  sound 
speed  profiles  for  750-Hz  runs.  Associated  environmental  data  are  given  in  Table  2.  A 
previously  conducted  acoustic  reflectivity  experiment  (4)  yielded  estimates  of  the  bottom 
sound  speed  and  density  in  the  area  as  1589  m/s  and  1.85  g/cm3,  respectively. 
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PROCESSING  TECHNIQUES 
Mode  Resolution  by  Time  of  Arrival 

If  a  duct  is  excited  with  a  pulse  of  duration  T0,  the  components  of  that  pulse  prop¬ 
agated  in  any  pair  of  normal  modes  of  the  duct,  having  group  velocities  Um  and  Un,  will 
be  completely  separated  in  time  at  range  (5) 

r>T0UmUnl(Um-Un).  (3) 

It  is  assumed  here  that  the  pulse  is  not  lengthened  significantly  by  dispersion.  Temporal 
resolution  was  used  to  examine  individual  modal  fields  whenever  group  velocities  differed 
sufficiently  to  produce  mode  resolution  at  a  range  where  the  available  signal-to-noise  ratio 
was  adequate.  When  complete  temporal  resolution  could  not  be  achieved,  a  spatial  pro¬ 
cessing  technique  was  used  to  resolve  the  modes. 


Mode  Resolution  by  Spatial  Filtering 

The  set  of  discrete  normal  modes  established  by  a  source  is  orthonormal.  That  is 

J1  Pi  «n  (?)  (?)  <*?  +  J°°  P2  un  (?)  «m  (?)  <*?  =  5m„>  (4) 

where  8mn  is  the  Kronecker  delta. 

The  orthogonality  permits  the  postulation  of  an  ideal  vertical  sensing  array  responding 
to  signals  propagating  in  any  one  selected  mode  and  rejecting  signals  propagated  in  any 
other  mode.  In  such  an  array  the  response  and  polarities  of  the  individual  elements  would 
be  adjusted  to  be  proportional  to  the  magnitudes  and  signs  at  corresponding  depths  of 
the  selected  mode  amplitude  function  (eigenfunctions).  An  ideal  array  would  be  capable 
of  monitoring  a  vertical  continuum  extending  from  the  water’s  surface  down  infinitely 
into  the  bottom.  The  sum  of  the  output  of  all  the  elements  of  such  rn  array  would  be 
zero  for  all  signals  except  those  propagating  in  the  selected  mode  as  indicated  in  Eq.  (4). 

In  practice  such  an  array  can  be  approximated  using  a  finite  number  of  sensors  (in  these 
experiments,  8  or  12)  spaced  throughout  the  water  column  and  adjusting  their  output 
gains  and  polarities  to  match  a  selected  mode  either  during  recording  (on  line)  or  pro¬ 
cessing  (playback).  The  princip  'l  limitations  on  the  effectiveness  of  this  array  are  imposed 
by  three  factors. 

The  first  is  interchannel  phase  distortion  resulting  from  anomalies  in  the  magnetic 
tape  recorder  (maybe  eliminated  by  on-line  processing)  and  from  bowing  of  the  line  array 
by  underwater  currents.  Such  distortion  degrades  the  coherent  summation  of  the  signals 
in  the  12  input  channels  resulting  in  a  reduced  filter  response  for  the  desired  mode.  Re¬ 
ductions  in  response  of  17%  to  22%  have  been  observed. 

The  second  arises  from  the  finite  number  of  discrete  samplings  of  the  signal  field. 
Although  the  modal  eigenfunctions,  in  theory,  are  orthogonal,  the  filter  may  not  be  com¬ 
pletely  effective  when  only  a  finite  sampling  of  the  field  is  available.  Sample  calculations 
for  the  array  used  in  these  experiments  indicated  that,  in  most  cases,  theoretical  rejection 
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of  at  least  10  dB  could  be  achieved.  Rejections  ranging  from  12  to  20  dB  have  been 
observed  experimentally  when  modes  were  temporally  resolved. 

The  third  is  induced  by  temporal  changes  in  the  mode  amplitude  functions.  This 
decreases  both  the  effective  rejection  of  unwanted  modes  and  the  signal  enhancement  of 
a  chosen  mode.  It  is  believed,  however,  that  this  was  not  as  impo.  -ant  a  factor  as  the 
first  two  because  there  were  indications  that  the  mode  amplitude  function  remained  rel¬ 
atively  stable  over  short  time  periods  (see  Temporal  Variations  section). 


Measurement  of  Vertical  Distributions  of  Pressure  Amplitude 

When  temporal  resolution  of  a  mode  was  achieved,  the  pressure  amplitude  of  the 
modal  field  was  measured  at  the  depth  of  each  of  the  12  hydrophones.  This  provided  an 
experimental  measure  of  the  distribution  of  pressure  amplitude  with  depth.  It  could  also 
provide  12  independent  measurements  of  the  mode  attenuation  coefficient  an,  since  all 
other  quantities  in  Eq.  (1)  are  known  or  can  be  calculated.  When  mode  resolution  re¬ 
quired  the  use  of  the  spatial  filter,  the  filter  output  was  used  to  calculate  an.  This  could 
be  done  since,  if  one  assumes  that  the  field  had  the  calculated  distribution  with  depth, 
the  filter  output  can  be  interpreted  either  in  terms  of  the  mean  signal  level  of  a  mode 
averaged  over  depth,  or  in  terms  of  the  modal  signal  level  at  any  given  depth. 
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RESULTS 

Comparison  of  Measured  and  Calculated  Vertical 
Distribution  of  Pressure  Amplitude 

A  normal  mode  computer  program  (3),  supplied  with  the  previously  stipulated  envi¬ 
ronmental  information,  calculated  theoretical  mode  amplitude  functions  (eigenfunctions), 
phase  and  group  velocities,  and  mode  wavenumbers. 

Experimental  vertical  distributions  of  pressure  amplitude  were  measured  for  five  runs 
in  which  modes  were  sufficiently  well  resolved.  The  theoretical  and  experimental  distri¬ 
butions  for  these  runs,  accompanied  by  the  sound  speed  profiles  used  to  calculate  the 
theoretical  mode  amplitude  functions,  are  shown  in  Figs.  6a  through  6e.  Only  the  first 
and  second  modes  had  sufficient  intensity  to  be  identified  and  measured.  Maximum  peak- 
to-peak  deflections  for  signals  propagated  in  those  modes  were  measured  from  each  of  the 
oscillographic  recorder  channels  for  ten  pulses  uniformly  distributed  throughout  a  run. 

The  average  of  the  ten  measurements  from  each  channel  was  plotted  at  the  normalized 
depth  of  the  hydrophone  associated  with  that  channel.  (A  phase  shift  of  approximately 
180°  was  observed  in  received  signals  whenever  the  pressure  amplitude  of  the  second  mode 
crossed  a  zero  value.)  The  resulting  experimental  mode  amplitude  values  and  polarities 
(dots)  were  then  displayed  with  the  calculated  mode  amplitude  functions  (solid  lines) 
after  an  amplitude  adjustment  using  a  scaling  constant  that  equated  their  maxima.  Three 
runs  were  from  the  March  1970  experiment  (750-Hz  source  frequency,  positive  gradient, 
and  nearly  isovelocity  sound  speed  profiles  —  Figs.  6a,  6b,  and  6c);  two  runs  were  taken 
from  the  July  1969  experiment  (400-Hz  source  frequency,  negative  gradient  sound  speed 
profiles  —  Figs.  6d  and  6e). 

Tbi  agreement  between  calculated  and  measured  mode  amplitude  functions  indicates 
that,  in  these  cases,  a  reasonably  accurate  prediction  of  mode  amplitude  distribution  with 
depth  was  possible.  The  slight  variations  in  the  measured  first-mode  amplitudes  from 
their  calculated  values  near  the  bottom  of  Fig.  6a  and  near  the  top  of  Fig.  6d  are  believed 
to  be  caused  by  the  arrival  of  the  third  mode.  The  group  velocity  of  the  third  mode  at 
the  source  center  frequencies  was  nearly,  but  not  exactly,  the  same  as  that  of  the  first 
mode  in  these  two  runs. 


Mode  Attenuation  Coefficients 

Mode  attenuation  coefficients  an  were  calculated  using  Eq,  (8)  and  measured  rn 
values  (see  Eq.  (7))  for  several  runs  for  400  Hz,  first  mode,  and  750  Hz,  first  and  second 
modes.  Intensities  were  calculated  for  each  pulse  from  the  adder  output  and  then  averaged 
to  obtain  a  time-averaged  value  of  intensity  for  a  run.  Table  3  shows  the  an  values  de¬ 
termined  from  the  time-averaged  intensity  and  velocimeter  cast  identifications  (these  ve- 
locimeter  cast  measurements  are  found  in  Ref.  2)  for  the  400-Hz  runs.  Since  intensity 
values  fluctuated  significantly  during  some  of  the  750-Hz  runs,  the  an  values  corresponding 
to  the  minimum  and  maximum  intensity  values  of  a  single  pulse  for  each  run  are  also 
displayed  with  the  time-averaged  <*n  values  in  Table  4.  Both  tables  contain  values  of  cx 
(1)  —  Cj  (0)  to  express  differences  in  the  water  layer’s  bottom  and  surface  sound  speeds 
for  the  prevailing  profiles.  (Since  all  profiles  were  essentially  monotonic  and  differences 
in  sound  speed,  if  any,  were  generally  localized  within  similar  depth  intervals,  these  values 
are  indicative  of  relative  strengths  of  the  sound  speed  gradients.)  The  tables  include  only 
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Fig.  6  —  Measured  and  calculated  vertical  distributions  of  pressure  amplitude  are  plotted  for  the  first 
and  second  modes  of  five  runs  in  which  these  modes  were  completely  separated  by  time  of  arrival. 
Also  included  are  plots  of  the  velocimeter  cast  measurements  associated  with  each  run. 
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Table  3 

Mode  Attenuation  Coefficients  at  400  Hz:  Time- 
Averaged  Values  for  Each  of  Several  Runs,  c  (1) 
—  c  (0)  Represents  Difference  Between  the  Water 
Layer’s  Bottom  and  Surface  Sound  Speeds. 


Mode  1 


Veloci- 

meter 

Cast* 

Ci  (D- 
Cl  (0) 
(m/s) 

“l 

(dB /km) 

11-3 

-17.3 

1.48 

6-1 

-16.0 

0.89 

1-6 

-15.8 

2.36 

2-1 

-15.8 

1.70 

6-5 

-15.7 

2.02 

11-5 

-15.0 

1.36 

12  2 

-14.5 

1.57 

9-6 

-13.6 

1.72 

*Velocimeter  casta  found  in  Ref.  2. 


Table  4 

Mode  Attenuation  Coefficients  at  750  Hz:  Minimum,  Maximum,  and  Time-Averaged 
Values  for  Each  of  Several  Runs.  Cj  (1)  —  c*  (0)  Represents  Difference  Between  the 
Water  Layer’s  Bottom  and  Surface  Sound  Speeds. 

Mode  1  Mode  2 


Veloci- 

meter 

Cast 

Cl  (1)- 
ci  (0) 
(m/s) 

aj  min 
(dB/ km) 

a  2  max 
(dB/ km) 

“1 

(dB/  km) 

a 

-3.0 

1.29 

1.64 

1.49 

c 

1.3 

-0.11 

0.37 

0.13 

d 

1.9 

-0.10 

0.25 

0.05 

e 

4.3 

0.07 

0.53 

0.28 

t 

4.4 

0.29 

0.52 

0.45 

g 

4.9 

0.14 

0.57 

0.35 

h 

6.1 

0.41 

1.00 

0.84 

i 

6.4 

-0.30 

-0.08 

-0.19 

j 

6.7 

0.03 

0.18 

0.12 

Veloci- 

meter 

Cast 

cj  (1>- 

c,  (0) 
(m/s) 

a2  min 
(dB/ km) 

a2  max 
(dr  km) 

“2 

(dB/ km) 

b 

-2.9 

0.53 

0.87 

0.64 

c 

1.3 

0.14 

0.72 

0.45 

d 

1.9 

0.37 

0.96 

0.77 

g 

4.9 

0.02 

0.91 

0.63 

g 

4.9 

0.45 

1.25 

0.92 

h 

6.1 

0.94 

1.41 

1.26 

j 

6.7 

0.62 

1.04 

0.85 

j 

6.7 

0.67 

0.90 

0.81 

. . r  imiiiir— Krai 
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those  data  obtained  from  runs  where  the  mode  eigenfunction  value  at  the  source’s  depth 
was  at  least  75%  of  its  maximum  value  in  the  water  column.  The  weaker  modal  signals 
resulting  from  less  mode  excitation  (see  Eq.  (1))  for  the  remaining  runs  were  judged  too 
small  to  be  confidently  distinguished  from  the  background  noise. 

Table  3  shows  generally  higher  first-mode  attenuation  values  for  400-Hz  than  for  the 
750-Hz  values  shown  in  Table  4.  Although  one  might,  in  general,  expect  a  decrease  in 
attenuation  for  a  lowering  of  frequency  (6),  it  is  believed  that  the  increased  bottom  inter¬ 
action  caused  by  the  very  sharp  negative  gradients  associated  only  with  the  400-Hz  runs 
were  sufficient  to  increase  their  attenuation  enough  to  obscure  the  effect  of  a  decreased 
source  frequency. 

Relations  among  mode  number,  sound  speed  gradients,  and  the  measured  attenuation 
maybe  observed  in  Table  4.  Primari'y,  an  increase  in  attenuation  accompanies  an  increase 
in  mode  number  for  all  positive  and  isovelocity  profiles.  The  negative-gradient  (profile  a), 
first-mode  run,  however,  had  the  highest  attenuation  in  Table  4.  Tne  first  mode  of  this 
run,  unlike  any  of  the  others  in  Table  4,  had  a  majority  of  its  energy  confined  to  propa¬ 
gate  in  the  slower  water  layer  near  the  bottom.  These  relations  parallel  those  that  might 
be  expected  if  only  that  portion  of  modal  attenuation  due  to  the  mode’s  interaction  with 
the  bottom  were  considered,  and  consequently  imply  a  strong  dependence  of  total  modal 
attenuation  on  modal  interaction  with  the  bottom. 

For  cases  in  which  modes  interact  relatively  weakly  with  the  bottom  (for  example, 
see  first-mode  attenuation  in  Table  4,  velocimeter  casts  c  through  ;),  other  attenuation 
mechanisms  might  be  expected  to  predominate.  Attenuation  for  confining  positive  gradi¬ 
ents  may  vary  because  of  interactions  with  the  surface,  such  as  scattering,  and  could  pos¬ 
sibly  depend  on  the  sea  state. 


Temporal  Variations 

Temporal  variations  of  the  modal  intensities  determined  from  the  mode  filter  output 
were  observed  in  the  750-Hz  runs  from  which  attenuation  coefficients  were  obtained  (see 
Table  4).  Plots  of  these  depth-averaged  intensities  vs  time  for  the  first  and  second  modes 
are  shown  in  Figs.  7a  and  7b,  respectively.  The  lower  case  letter  by  each  run  designates 
the  velocimeter  cast  associated  with  that  run  (see  Table  1).  These  temporal  variations 
could  be  caused  by  changes  in  the  pressure  amplitude  distribution  or  by  fluctuations  of 
mode  attenuation  or  excitation.  A  test  was  conducted  in  an  attempt  to  clarify  the  source 
of  the  temporal  variations. 


In  this  test,  temporal  fluctuations  of  the  amplitude  of  modal  arrivals  at  individual 
hydrophones  were  compared  with  the  temporal  fluctuations  of  the  filter  output.  The 
maximum  amplitude  of  each  signal  propagated  in  a  particular  mode  was  measured  for  a 
sequence  of  approximately  60  pulses,  both  at  the  output  of  an  individual  hydrophone 
and  at  the  output  of  the  mode  filter  matched  to  that  mode.  A  cross-correlation  coeffi¬ 
cient  was  then  calculated  for  those  sequences  using 
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where 

Rn  =  nth  mode’s  cross-correlation  coefficient 
I  =  total  number  of  pulses  in  a  run 

and  Vj  =  outputs  of  an  individual  hydrophone  and  the  mode  filter,  respectively, 
for  the  ith  pulse 

x  and  y  =  arithmetic  means  of  those  outputs  over  /  pulses. 

When  this  was  done  for  several  hydrophones  in  each  of  two  runs  the  correlation  coeffi¬ 
cients  for  the  first  mode  ranged  from  0.72  to  0.92.  For  the  second  mode  they  ranged 
from  0.56  to  0.92.  The  correlations  obtained  tend  to  indicate  that  at  least  the  major 
features  of  the  temporal  fluctuations  observed  at  discrete  depths  coincided  with  those  of 
the  field  averaged  over  depth.  Thus,  it  appears  that  the  major  cause  of  the  observed  fluc¬ 
tuations  is  associated  with  changes  in  total  field  strength  rather  than  with  changing  dis¬ 
tributions.  Further  evidence  of  this  can  be  observed  in  Figs.  6a  through  6e  where  the 
time-averaged  measured  distributions  agree  well  with  the  calculated  distribution. 

On  the  other  hand,  the  very  large  (up  to  9  dB  at  10  km  in  7  min)  changes  in  signal 
level  shown  in  Figs.  7a  and  7b  are  difficult  to  account  for  through  variable  attenuation 
mechanisms.  However,  one  might  expect  the  vertical  pressure  amplitude  distribution  of 
the  measured  modes  to  be  more  strongly  dependent  on  environmental  conditions  in  the 
vicinity  of  the  receiving  hydrophones.  If  this  were  the  case,  environmental  conditions 
more  distant  from  those  hydrophones  might  affect  the  average  modal  signal  amplitudes 
without  significantly  altering  the  distributions. 


CONCLUSIONS 

Comparisons  of  measured  and  calculated  vertical  distributions  of  pressure  amplitude 
for  the  first  two  modes  of  400-  and  750-Hz  signals  demonstrated  that  such  distributions 
could  be  measured  and  predicted  in  at  least  one  simple  experimental  geometry  for  various 
sound  speed  profiles  including  positive,  zero,  and  negative  gradients.  It  remains  to  extend 
these  capabilities  to  areas  of  different  bottom  types  and  more  complicated  bottom  contours. 

The  prediction  of  the  pressure  amplitude  distribution  with  depth  proved  sufficiently 
accurate  to  encourage  the  development  of  a  spatial  mode  filtering  technique  based  on  the 
inherent  orthogonality  of  normal  modes.  Such  a  filtering  technique  was  used  to  determine 
Individual  mode  intensities  averaged  over  depth. 

Mode  intensity  variations  of  up  to  9  dB  at  a  range  of  10  km  were  observed  within 
7-min  runs.  Specific  causes  of  these  variations  are  not  known.  Various  cross-correlation 
tests  tended  to  indicate,  however,  that  the  vertical  distributions  of  pressure  amplitude 
remained  relatively  stable.  Therefore,  fluctuations  in  attenuation  mechanisms  that  affect 
the  total  signal  field  uniformly  with  depth  at  the  receiving  line  of  hydrophones  are  sus¬ 
pected  to  be  the  predominant  sources  of  the  observed  mode  intensity  variations. 

The  average  attenuation  values  of  the  first  two  modes  for  several  runs  were  determined 
from  the  mode  intensity  measurements.  Relations  among  mode  numbers,  sound  speed 
gradients,  and  attenuation  values  showed  an  observable  dependence  of  the  modal  attenua¬ 
tion  on  the  degree  to  which  the  mode  interacts  with  the  bottom.  (Attenuation  increasing 
with  increasing  bottom  interaction.)  For  situations  where  positive  gradients  are  suffi¬ 
ciently  sharp  to  confine  modes  to  water  layers  near  the  top,  attenuation  mechanisms  other 
than  those  associated  with  the  bottom  may  become  prevalent. 
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